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We introduce a necessary and sufficient criterion for the non-Markovianity of Gaussian quantum dynamical 
maps based on the violation of divisibility. The criterion is derived by defining a general vectorial representa¬ 
tion of the covariance matrix which is then exploited to determine the condition for the complete positivity of 
partial maps associated to arbitrary time intervals. Such construction does not rely on the Choi-Jamiolkowski 
representation and does not require optimization over states. 
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In recent years much effort has been devoted to the charac¬ 
terization and quantification of non-Markovianity in the evo¬ 
lution of open quantum systems (see e.g. Ref. for a re¬ 
cent review). Non-Markovian quantum evolutions may typ¬ 
ically arise in the presence of structured environments, such 
as in quantum biological systems Hi] and in squeezed baths 
of light with finite bandwidth js]]. Moreover, recent studies 
suggest that properly engineered non-Markovian channels can 
improve the efficiency of quantum technology protocols due 
to the backflow of information from the environment to the 
system sni. Establishing whether noisy quantum evolu¬ 
tions are non-Markovian and therefore preserve some mem¬ 
ory on the story of the system is of capital importance in the 
field of quantum cryptography ifT^ . 

Various approaches to the characterization and quantifica¬ 
tion of quantum non-Markovianity have been introduced in 
recent years UBS . Most of them are witnesses, and thus 
rely on sufficient, but not necessary, conditions H. They usu¬ 
ally are based on the non-monotonic behavior of certain quan¬ 
tities in the presence of memory effects iflbl - l^ . Moreover, 
most of them rely on optimization over states. 

Proper measures of non-Markovianity have also been in¬ 
troduced for finite-dimensional systems. These measures in¬ 
clude the amount of isotropic noise necessary to make the dy¬ 
namics completely positive in every arbitrary short interval of 
time 123] and the negativity of the decay rates appearing in 
the generators of the time evolution, once the associated mas¬ 
ter equation is expressed in canonical form ll^ . 

A further necessary and sufficient criterion has been ob¬ 
tained by Rivas, Huelga, and Plenio (RHP) by considering 
the violation of the divisibility property, which expresses the 
possibility of decomposing the evolution on a generic time 
interval into two successive, independent completely posi¬ 
tive maps. Non-Markovianity is then characterized by the 
extent that the intermediate map violates complete positivity 
(CP) d. These three necessary and sufficient criteria for 


‘Corresponding author: illuminati@sa.infn.it 


finite-dimensional systems have been shown to be completely 
equivalent JUt]; moreover, the criterion based on the isotropic 
noise and the RHP criterion rely on positivity of the Choi- 
Jamiolkowski states corresponding to the channels ll^l2^ . 

Addressing the general characterization and quantification 
of non-Markovianity in the infinite-dimensional case is an 
important open question, given the great importance of, for 
instance, Gaussian states and Gaussian channels in quantum 
optics, quantum information, and quantum technologies. In¬ 
spired by the RHP approach, in the present work we introduce 
a necessary and sufficient criterion of non-Markovianity for 
Gaussian evolutions. It is based on the violation of the divis¬ 
ibility property directly at the level of the matrices defining 
the channels, exploiting a powerful vectorial representation 
for them and for the covariance matrix. As such, the criterion 
does not require the use of the Choi-Jamiolkowski isomor¬ 
phism between states and channels ES. Moreover, it does 
not require optimization over states, a chal leng ing task both 
for finite- and infinite-dimensional systems 11^ . 

Based on such criterion, we introduce the corresponding 
measure of non-Markovianity for Gaussian channels and il¬ 
lustrate it for some paradigmatic examples. For the specific 
channels considered, violation of divisibility turns out to be 
equivalent to the negativity of the decoherence rates appear¬ 
ing in the canonical form of the master equation. 

Given a generic input state, its time evolution in a Gaussian 
channel is defined according to the following transformation 
on the input covariance matrix ct(0): 

<j{t)=X{t)cj{Q)X^ {t) + Y{t) , (1) 

where (X, Y) are 2N x 2N real matrices; moreover Y is sym¬ 
metric. It is possible to show that the CP requirement imposes 
the condition : 

Y{t)-Ui+'-x{t)nx^ {t)>Q , ( 2 ) 

where Ul is the symplectic matrix, and the symbol j denotes 
matrix transposition. 

Gaussian channels enjoy a semigroup structure. Given two 
such channels corresponding, according to Eq. ([ill, to the 
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pairs {Xi,Yi) and th e resulting composed channel 

is characterized as follows 

(Xi,ri) • {X2,Y2) = iX^X2,XiY2Xj + Yi) . (3) 

Consider now a quantum evolution from time tg to t 2 de¬ 
scribed, in general, by some family of trace-preserving linear 
maps {£{t 2 , ti), t 2 > ti > to}. The map is said to be divisi¬ 
ble, or Markovian, if, for every t 2 and ti it holds that 

£it2,to) = £it2,ti)£(tl,to), t2>ti>to, (4) 

and the map f (^ 2 , f i) is CP. An evolution is non-Markovian if 
it violates the divisibility property, Eq. (|4]i. 

For Gaussian channels, we can reformulate the general di¬ 
visibility condition, Eq. in the following form. Let us first 
introduce an auxiliary vectorial notation. In this notation the 
elements of the covariance matrix cr form a vector according 
to a lexicographical ordering, i.e. {a{t))i^ = a{t)ij (where 
k = N{i — 1) + j and i,j = 1,..., 2N), which in the more fa¬ 
miliar Dirac notation can be written as {k\(j{t)) = {ij\a{t)) = 
{i\a{t)\j). For convenience, we also add to a an auxil¬ 
iary vector entry of value 1. In this notation one can obtain 
the following representation; [X(t)a(0)X'^(t)]ij =[(2f(f)0 
X(fpCT(0)]fc = [<I>(f)CT(0)]fc, where <I>(f) = 2f(f) (g) X{t). 
In Dirac notation, one has: 

(i|X(f)cr(0)XT(f) \j) = {i\X{t)\n) (n|a(0) |m) {m\X'^{t) \j) 

n,m 

=y^{ij\X{t)^X{t)\nm){nm\(7(0)) = {ij\X(t)'^X{t)\(7{0)), 

n^m 

where 1^) (^1 '^^e identity resolution in some basis. One 

can now reexpress Eq. ([T]i in terms of a vector by matrix mul¬ 
tiplication; 

(,-(„) ?(„)(,-( 0 ))_ ,,, 

where 0 = (0,..., O)'*' is the 2A^-dimensional null vector and 
Y (t) is the vectorial form of the matrix Y(t). Vectorization is 
an isomorphism, thus reversible: de-vectorizing Eq. 0 yields 
exactly the standard representation, Eq. O. 

In vectorial notation, the channel composition law, Eq. 0, 
is reexpressed by the following matrix multiplication form; 

( 0" 1 j I OT 1 j = ( OT 1 j ■ ® 

Setting for ease of notation to = 0, G = t, and t 2 = t + e 
for any instance of time t and e, by the continuity of time the 
dynamics can be split as [0, f -f e] = [0, f] U [f, f -I- e], and one 
can obtain the vectorial expression for the intermediate map 
in the interval [t, t + e]. Let us comment on the invertibility 
of the X{t,0) matrix. Examples of Gaussian channels char¬ 
acterized by a non invertible X matrix can be found based 
on classification of one-mode Gaussian channels provided in 
Ref. iH. Up to Gaussian unitary equivalence, channels for 


which X is non invertible include the completely depolaris¬ 
ing channel which projects every input state on a thermal state, 
and channels which transform the canonical quadrature Q and 
P as: P ^ p, Q ^ Q + q, where p and q are thermal states. 
However, non invertible cases do not impose any restriction 
on our procedure, because one can always introduce the ma¬ 
trix lr] + X{t, 0), determine its invers^ and evaluate the limit 
p ^ 0, which is always non-singular IUI gIII^ . 

De-vectorizing the intermediate Gaussian map, we obtain 
its complete expression in terms of the X and Y matrices; 

X{t+ e, t)=X(t+e, Q)X-\t, 0), 

(7) 

Y{t+e, t)=Y{t+e, 0)-X{t+€, t)Y{t, 0)X'^{t+e, t). 

Since the condition of divisibility is equivalent to CP of the 
intermediate map, Eq. 0, fromEqs. 0, 0, and 0, the con¬ 
dition for non-Markovianity at any given time t reads: 

Y{t + e,t)-Ul + ^X {t + e, t) UXT {t + e,t)<0. (8) 

Given that Eq. 0 is necessary and sufficient for CP it 
follows that Eq. 0 is a necessary and sufficient criterion for 
the non-Markovianity of Gaussian channels. 

The necessary and sufficient condition, Eq. 0, allows to 
introduce a proper measure of non-Markovianity for Gaussian 
channels by quantifying the extent by which the intermediate 
dynamics fails to be CP. This corresponds to the quantification 
of the negative part of the spectrum of the symmetric matrix 
appearing on the l.h.s. of Ineq. 0. Denoting the set of eigen¬ 
values by {vk{t + e, t)}k=i,...,2N, the following functions 

fk{t) = \ lim [\vk{t + e,t)\- Vk{t + e,t)] (9) 

Z e—>-0+ 

quantify the negative contribution at time t given by the fcth 
eigenvalue. Therefore punctual non-Markovianity, quantified 
by the negative part of the spectrum at a given time t, reads 

2N 

F{t) = J2hit) . ( 10 ) 

k=l 

Since F{t) > 0 if and only if the evolution is non-Markovian, 
and F{t) =0 otherwise, total non-Markovianity on a generic 
time interval I is 



It is important to note that, when the dynamics is described 
by means of a master equation, the expression of the matrices 
{X,Y} that define the channel are obtained directly, in the 
phase space formalism, through Eq. O, from the expression 
of the characteristic function of the evolved Gaussian state. 

The above divisibility-based necessary and sufficient crite¬ 
rion is completely general: it holds for any Gaussian map, 
independently of the existence of a generator. On the other 
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hand, as already mentioned. Hall, Cresser, Li, and Anders- 
son have recently shown that in the finite-dimensional case, 
for which, at variance with the infinite-dimensional case, all 
processes always admit a generator, the necessary and suffi¬ 
cient criterion for non-Markovianity based on divisibility is 
equivalent to the criterion based on the negativity of the de¬ 
coherence rates appearing in the canonical form of the master 
equation It is then tempting to conjecture that this equiv¬ 
alence holds also in the infinite-dimensional case for channels 
that admit a generator, we discuss two paradigmatic cases that 
admit a representation in terms of canonical master equations, 
pure damping and quantum Brownian motion, and show that 
in such instances the equivalence indeed holds. 

The simplest example we can begin with is the Lindblad- 
type master equation describing the damping process for a 
single field mode with a single decay rate: 


dp jt) 
dt 


a 7(f) 





( 12 ) 


where a <C 1 is the coupling constant and 7 (f) is the damping 
rate. The evolution of a generic Gaussian state in this Gaus¬ 
sian channel is described by the corresponding evolution of 
the displacement and covariance matrices. From the latter, via 
Eq. O, one obtains the X and Y matrices: 


. X r(t) 

X{t,0) = e“ —1 , 

(13) 

y(f,0)= [i-e-r«] i , 

(14) 


where r(t) = 2a J^j(s)ds. Eqs. (fT3]) . (fT4li allow to ob¬ 
tain, through Eqs. (|7]), the matrix appearing in the l.h.s. of 
the CP condition, Ineq. It is straightforward to ver¬ 

ify that the eigenvalues of this matrix are negative when 
exp(— r(t-|-e, t)) < 1 , where r(f-|-e,f) = r(t-|-e,0)—r(f,0). 
Moreover, to first order in e, we have r(t + e,t) ss 2j{t)e. As 
a consequence, the evolution is non-Markovian if and only if 
7 (t) < 0 , showing that in this case violation of divisibility is 
indeed equivalent to negativity of the decoherence rate. The 
corresponding measure, through Eqs. (|9l) and (flOl) reads: 


=—aJ 7 (f)(if, (15) 

where /' are the sub-intervals of I in which 7 (f) < 0. 

We next consider Quantum Brownian Motion in the weak 
coupling limit and under the secular approximation. It is de¬ 
scribed in the interaction picture by the following Lindblad- 
type master equation (see Ref. S' and references therein): 

dp (f) A (f) -I- 7 (f) f ^ ,, 

^ [2apat - {a^a,p}] + 

7(i) [ 20 ^^ - {aa^jo}] ■ (16) 

The coefficients 7 (f) and A(f), in general time-dependent, are 
respectively the damping and diffusion coefficient, whose ex¬ 
plicit expressions are obtained once one selects the explicit 
form for the spectral density of the bath. The general solution 


of Eq. (fThl) allows to obtain the evolution of the displacement 
and covariance matrices for any input Gaussian state 1^ . The 
corresponding X and Y matrices read: 

r(t) 

X{t,0) = e- — R{t) , (17) 

r(f,0) =e"’"(*^A(f)l , (18) 


where r(f) = 2 f*j(s)ds, A(f) = f* e^^^^A(s)ds, R{t) is 
the rotation matrix by the angle tuot, and ujq is the system’s 
characteristic frequency. These expressions and Eqs. © de¬ 
termine the eigenvalues of the matrix in the l.h.s. of Ineq. (| 8 ]i: 


vi (f -h e,f) 
i^2{t + e,f) 


g-r(t+e,t) 4. 2 A (f -h e, f) - 1 , 

1 - +2A{t + e, t) e-r(‘+«.o) ^ 

(19) 


where r(f -|- e, f) = r(f -|- e, 0) — r(f, 0) and A(f -I- e, f) = 
A(f-|-e, 0) —A(f, 0). To first order in e, we have: ~ 

1 — 27 (f)e and A(f -I- e,f) ~ A(f)e. Then, condition 
Eq. ([ 8 ]l on the eigenvalues, i.e. the violation of the divisibil¬ 
ity condition, implies A(f) < | 7 (f)|. This is again equiv¬ 
alent to negativity of the decoherence rates [A(f) -|- 7(f)]/2 
and [A(f) — 7(f)]/2 appearing in Eq. (fT 6 l l. Einally, exploiting 
Eqs. ( fT^ and I© we obtain the following expression for the 
punctual measure of non-Markovianity: 

E’(f) = i[|A(f)-7(f)| + |A(f)+7(f)|]-A(f) . (20) 

In order to investigate explicitly the behavior of non Marko- 
vianity in the Quantum Brownian Motion, we need to spec¬ 
ify the spectral density to obtain explicit expressions of the 
damping and diffusion coefficients 7 (f) and A{t). Consider¬ 
ing the rather typical case of an Ohmic bath with an exponen¬ 
tial cut-off uJc, the parameters that govern the dynamics are 
the temperature T and the ratio between the cut-off frequency 
of the bath and the characteristic frequency of the system 
X = iOcjiOQ. It is expected that, in the regime x <C 1 the dy¬ 
namics should be non-Markovian, while Markovianity should 
be recovered for x ^ 1 1^ . It is also convenient to ex¬ 
press the evolution in terms of the dimensionless reduced time 
T = uJct. Moreover, explicit analytic expressions of the diffu¬ 
sion coefficient A(t) can be obtained quite straightforwardly 
in the high- and low-temperature regimes ll^ . Considering 
first the asymptotic values of the damping and diffusion co¬ 
efficients in the large-time limit, t —> 00 , both in the high- 
and low-temperature regimes, it is straightforward to verify 
that the asymptotic punctual non-Markovianity E’(c)o) = 0: at 
large times Markovianity is always recovered, independently 
of the values of the parameters that govern the dynamics. 

Considering now generic times, in Eig. ([Tlr) we report the 
behavior of the punctual non-Markovianity F, Eq. ( l20b . as a 
function of the reduced time t at fixed values of the parame¬ 
ter X = oJc/ojQ in the high-temperature limit. In this regime 
A(r) ^ 7 (r) lO . and the non-Markovianity of the dy¬ 
namics depends essentially only on the diffusion coefficient: 
F{t) ~ |A(r)| — A(r). Hence, the time interval for which 
















4 




FIG. 1: (color online) Non-Markovianity F{t). a): High- 

temperature limit for a; = 0.1 (blue full line), a: = 0.2 (red dot- 
dashed line), and x = 0.3 (green dashed line). Inset: A(r) for 
a; = 0.1 (blue full line), x = 0.2 (red dot-dashed line) and x = 0.3 
(green dashed line). 6): Low-temperature limit for x = 0.2 (blue full 
line), a; = 1.0 (red dot-dashed line) and x = 2.0 (green dashed line). 
Inset: A(r) (red dot-dashed line) and 7 (r) (blue line) for a; = 0.1. 


the evolution is non-Markovian, F{t) > 0, corresponds to 
the negativity of the decoherence rate, A(r) < 0. Non- 
Markovianity is strongest in the regime a; <C 1, corresponding 
to the characteristic time of the bath being much larger than 
the characteristic time of the system. When x increases, the 
negative part of the oscillations and F{t) quickly vanishes, 
and one recovers the Markovian regime. 

In the low-temperature regime, see Fig. ([TI)), the diffu¬ 
sion and damping coefficients are comparable, and the non- 
Markovianity F{t) is given by the full expression, Eq. ( |20] |. 
In this situation, a non-Markovian regime is observed also for 
A(r) > 0, provided A(r) < 7 (r), and even if the character¬ 
istic times of the bath start to be comparable or smaller than 
the characteristic times of the system, x > 1. 

In these examples, the criterion based on the X and Y ma¬ 
trices defining a Gaussian channel turns out to correspond to 
the negativity of the decoherence rates. On the other hand, it 
should be stressed that the criterion is much more general and 
applies to any Gaussian evolution, including those that do not 
admit a generator and hence cannot be described in terms of 
master equations. Finally, it always allows, at least in prin¬ 
ciple the experimental verification of the Markovianity of the 
evolution. A further advantage is that such verification does 
not require optimization over the set of input states, since it is 


based directly on the characteristic matrices that define intrin¬ 
sically the dynamical map. When the generator exists, so that 
the dynamical map can be associated to a master equation, the 
general quantifiers of non-Markovianity, Eqs. (|9l) and (fTTl i. re¬ 
duce to simple functions of the decoherence rates, Eq^dni 
and ( l20l i. which can be reconstructed experimentally O^ . 

Summarizing, checking for CP of the partial map in the 
finite-dimensional case consists in verifying the positivity of 
the corresponding Choi-Jamiolkowski state. The RHP mea¬ 
sure of non-Markovianity is then defined in terms of the neg¬ 
ative part of the spectrum of such state din. One might 
try to enforce the same criterion, in complete analogy with 
the finite-dimensional case, by checking CP of the interme¬ 
diate Gaussian map by checking the positivity of the corre¬ 
sponding Choi-Jamiolkowski state. Indeed, for single-mode 
CP Gaussian maps, it has been shown that a Kraus decompo¬ 
sition can always be found, so that one can always construct 
the corresponding Choi-Jamiolkowski state 1^ . Unfor¬ 
tunately, the Kraus representation does not exists for non CP 
Gaussian maps and the Choi-Jamiolkowski states correspond¬ 
ing to these maps are to date not characterized. Therefore, 
checking violation of CP for Gaussian maps using the Choi- 
Jamiolkowski isomorphism is currently impossible. 

We have succeeded in circumventing this stumbling block 
by expressing the condition of CP of general Gaussian maps, 
Ineq. (I2]i, directly in terms of the X and Y matrices govern¬ 
ing the evolution of the covariance matrix of input Gaussian 
states. Concatenation of the maps in the finite-dimensional 
case is expressed straightforwardly by matrix multiplication, 
thanks to the superoperator representation of the quantum 
channel. Immediate generalization of this method to general 
Gaussian maps is not possible. Instead, we succeed in defin¬ 
ing Gaussian state vectorization by a suitable vectorization of 
the covariance matrix. This allows to introduce a matrix rep¬ 
resentation also for Gaussian maps and generalize the RHP 
method of characterizing non-Markovianity. 

The specific form of the vectorization procedure that we 
have introduced is suitably defined in such a way to preserve 
the fundamental semigroup property of Gaussian channels, al¬ 
lowing to investigate any Gaussian map in a compact and ele¬ 
gant form. 

Since our approach does not require optimization on the 
set of input states, it can be especially helpful when consider¬ 
ing multi-mode Gaussian channels, allowing in principle for a 
systematic study of the interplay between non-Markovianity, 
entanglement, and coherence. Furthermore, it can open the 
way to the characterization of Gaussian quantum metrol¬ 
ogy in non-Markovian environments, extending to Gaus¬ 
sian states of continuous-variable systems the investigations 
pioneered in Refs. 
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